Polaronic signatures and spectral properties 
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We explore the consequences of electron-phonon (e-ph) coupling in graphene antidot lattices 
(graphene nanomeshes) , i.e., triangular superlattices of circular holes (antidots) in a graphene sheet. 
They display a direct band gap whose magnitude can be controlled via the antidot size and density. 
The relevant coupling mechanism in these semiconducting counterparts of graphene is the modula- 
tion of the nearest-neighbor electronic hopping integrals due to lattice distortions (Peierls-type e-ph 
coupling). We compute the full momentum dependence of the e-ph vertex functions for a number 
of representative antidot lattices. Based on the latter, we discuss the origins of the previously found 
large conduction-band quasiparticle spectral weight due to e-ph coupling. In addition, we study 
the nonzero-momentum quasiparticle properties with the aid of the self-consistent Born approxima- 
tion, yielding results that can be compared with future angle-resolved photoemission spectroscopy 
measurements. Our principal finding is a significant e-ph mass enhancement, an indication of pola- 
ronic behavior. This can be ascribed to the peculiar momentum dependence of the e-ph interaction 
in these narrow-band systems, which favors small phonon momentum scattering. We also discuss 
implications of our study for recently fabricated large-period graphene antidot lattices. 

PACS numbers: 63.20.kd, 63.22.-m, 65.80.Ck, 73.21. Cd 
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I. INTRODUCTION 

The discovery of freestanding graphene has spawned 
huge interest in the electronic and transport proper- 
ties of this material.^ In particular, a great deal of re- 
search effort is presently being dedicated to graphene- 
based superlattices^^S Among them, a class of superhon- 
eycomb systems^ - graphene antidot lattices {graphene 
nanomeshe^) - was introduced^ by analogy with the 
concept of antidot lattices defined atop a two-dimensional 
electron gas in a semiconductor heterostructure.^ The 
electronic fi2r22 optical, and magnetioi^ properties of 
these superlattices result from a subtle interplay of the 
intrinsic features of graphene and a lattice periodicity 
imposed by holes (antidots) in a graphene sheet. 

Aside from some fundamental aspectSfi^ the main in- 
centive behind the current graphene-related research ac- 
tivity stems from the prospect of a carbon-based elec- 
tronics.^^ Graphene displays exceptional properties, such 
as room-temperature ballistic transport^ on a submi- 
cron scale with intrinsic mobility as high as 200, 000 
cm^/Vs and the possibility of heavy doping without alter- 
ing significantly the charge-carrier mobility. Yet, intrinsic 
semimetallic graphene is of limited utility for fabricating 
electronic devices!^ because the transmission probabil- 
ity of Dirac electrons across a potential barrier is always 
unity - independent of the height and width of the barrier 
- a manifestation of Klein tunneling.^* As a result, the 
conductivity cannot be controlled through a gate volt- 
age, which, however, is a prerequisite for operation of 
a conventional field-effect transistor where switching is 
achieved by gate-voltage induced depletion. Thus it is 
crucial to open a band gap in graphene. A possible ap- 
proach entails processing graphene sheets into nanorib- 



bons, which show a confinement-induced gap that scales 
inversely with their width. To get a gap sufficient for 
room-temperature transistor operation (Eg ^ 0.5 eV) 
one needs ribbons with a width of only few nanometers. 

Graphene antidot lattices have a direct electronic band 
gap stemming from the quantum confinement brought 
about by the periodic potential of a regular arrangement 
of antidots. The magnitude of the gap follows a generic 
scaling relatior^i^ - a linear increase with the product of 
the antidot size and density - and can thus be exter- 
nally controlled. It is this tunability of the band gap - 
not available when a gap is induced by growing graphene 
epitaxially on substrates such as SiC - that makes these 
superlattices interesting from both the fundamental and 
practical viewpoints. They can be thought of as highly 
interconnected networks of nanoribbons, and can in prin- 
ciple be fabricated on suspended graphene monolayers 
via electron-beam lithography, whose current resolution 
limit is as low as 10 nm. Yet, being serial in nature, 
this method is not applicable to large-area patterning 
of graphene. An alternative fabrication approach, block 
copolymer lithography, has recently been utilized to fab- 
ricate very uniform, free-standing triangular graphene 
antidot lattices.— "iS Superior room-temperature electri- 
cal properties of the fabricated devices - compared to 
graphene nanoribbons - have also been demonstrated 

A thorough understanding of e-ph scattering is of ut- 
most importance because the latter determines the ul- 
timate performance limit of any electronic device: in 
the high-bias regime of transport the e-ph scattering in- 
creases the differential resistance. In the present work, we 
address the e-ph coupling^*' in graphene antidot lattices 
and explore some of the ensuing physical consequences. 
The major e-ph-interaction mechanism in this system, 
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as well as in 7r-electron systems in general)2i"— is the 
modulation of electronic hopping integrals due to lattice 
distortions (Peierls-type e-ph coupling^i). 

It is by now widely accepted that e-ph coupling 
in graphene is comparatively weaki ^^'^^ for instance, 
the angle-resolved photoemission spectroscopy (ARPES) 
data on inelastic carrier lifetime^! were consistently ex- 
plained without even invoking phonon-related effects.— 
Here we demonstrate that in graphene antidot lattices, 
which have completely different band structure than their 
"parent material", e-ph coupling is very significant, i.e., 
properties of these systems undergo a strong phonon- 
induced renormalization. We ascribe this not only to 
their narrow-band character and low dimensionality, but 
also to the peculiar momentum-dependence of the e-ph 
interaction, favoring small phonon momentum scattering. 

Unlike in our recent work,— which was solely con- 
cerned with the zero-momentum conduction-band quasi- 
particle spectral weight for different graphene antidot 
lattices, in the present work we put an additional em- 
phasis on the directional dependence of the e-ph mass 
enhancement. This is particularly appropriate because 
of the anisotropic (in momentum space) character of the 
relevant e-ph coupling. We also describe the nonzero- 
momentum spectral properties due to the e-ph coupling 
at the level of the self-consistent Born approximation. In 
this way, we obtain results that can be compared with 
future ARPES measurements. 

The outline of the paper is as follows. In Sec. |ll] we in- 
troduce the system under consideration, along with the 
notation and conventions to be used throughout. In sub- 
sequent Sees. Illll and llVl we briefly recapitulate the main 
features of the band structures and phonon spectra of 
graphene antidot lattices, respectively. Sec.|V]is set aside 
for a detailed account of the microscopic mechanism of e- 
ph coupling in the system, with emphasis on its resulting 
momentum dependence. The following Sec. lVIl is devoted 
to the effect that this momentum-dependence has on the 
phonon-induced renormalization of electronic properties. 
In Sec. IVIII we discuss some implications of our study for 
optical absorption experiments and charge transport, as 
well as possible future work. We conclude, with a brief 
summary of the paper, in Sec. I Villi Some mathematical 
details are deferred to Appendices E] and [BJ 



II. NOTATION AND CONVENTIONS 

Triangular graphene antidot lattices with circular anti- 
dots are illustrated in Fig. [TJ Their geometry is uniquely 
specified by two dimensionless parameters: the side 
length of their hexagonal unit cells (L) and the antidot 
radius (i?), both expressed in units of the graphene lat- 
tice constant a ~ 2.46 A. (Note that a = accVi, where 
Qcc = 1-42 A is the distance between nearest- neighbor 
carbon atoms.) While L is invariably an integer, and 
defines the superlattice period (center-to-center 
distance between nearest- neighbor antidots), R can also 




FIG. 1: (Color online) A finite segment of the graphene an- 
tidot lattice {L, R} with circular antidots. Here ai and sl2 
are the superlattice basis vectors; L and R are, respectively, 
the side length of the hexagonal unit cell and the antidot ra- 
dius, both expressed in units of the graphene lattice constant 
a ~ 2.46 A. The vectors di, 62, and 63 specify positions of the 
nearest neighbors of a carbon atom on sublattice A. 



take non-integer values. Therefore, the notation {L, R} 
will hereafter be used to specify different antidot lattices. 
If a carbon atom (henceforth abbreviated as C atom) on 
sublattice A is taken to be the origin, the adjacent C 
atoms are determined by the vectors Si, 62, and ^3 (see 
Fig.[T|). (For an origin at sublattice B, the corresponding 
vectors would be —Si,— 82,— S3.) 

In what follows, the vectors R will designate the 
unit cells {N of them) of an antidot lattice, while dm 
(to = I, . . . , A^at) will specify the positions of the C atoms 
within a unit cell. Thus the vectors R + dm uniquely 
specify the positions of all C atoms. 

III. BAND STRUCTURE OF GRAPHENE 
ANTIDOT LATTICES 

Much like some other types of graphene-based super- 
lattices,^ typical antidot lattices contain extremely large 
numbers of C atoms in their unit cells; in particular, here 
we consider superlattices with A''at ~ 300 — 1600 atoms 
per unit cell. Such a large system size prohibits the use of 
methods based on the density functional theory (DFT) 
for the band-structure and phonon calculations. Instead, 
we model the band structure using a nearest-neighbor 
TT-orbital tight-binding model Hamiltonian 

-^e = J2 (ak+d„+5«R+d„ + H.c. ) , (1) 

R,, m,5 

where aR+d„ ('^R+d ) destroys (creates) an electron in 
the 2pz orbital (p{r — R — dm) of the C atom located at 
the position R + dm , S stands for the nearest neighbors 
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of that C atom, and t « 2.8 eV is the nearest-neighbor 
hopping integral. The tight-binding method is known 
to be capable of reproducing very accurately the low- 
energy part of the DFT band structure of graphene^i Its 
accuracy in the case of graphene antidot lattices, i.e., its 
good agreement with the DFT results for lattices with 
very small unit cells, has recently been demonstrated.^^ 
The Bloch wave functions corresponding to the energy 
eigenvalues en(k), where n is the band index and k is a 
quasimomentum in the Brillouin zone (BZ), are given by 

V-nkW^^C'V^kW, (2) 

m 

where (/)mk(r) is the discrete Fourier transform of fir — 
R— d„). Due to the sixfold-rotational point-group sym- 
metry of the system, the coefficients C^'*' obey the con- 
ditions 

/On,k „ie(k)/^",k' fn\ 

'^m — ^ '^m' ' \'^) 

where 0(k) is an irrelevant k-dependent phase and k' 
(m') is obtained by rotating the vector k (atom m) 
through an angle of vr/S. 

Unlike graphene itself, graphene antidot lattices show a 
band structure characteristic of direct band gap semicon- 
ductors. In particular, the lattices with circular antidots, 
that we are concerned with here, fall into the class of su- 
perhoneycomb systems (with a unit cell that must have 
at least the six-fold rotational symmetry) whose possible 
band structures were classified by Shima and Aokiji 

In what follows, we study the {L, i? = 5} (9 < L < 19) 
and {L, i? = 7} (12 < L < 20) families of graphene 
antidot lattices. We find that these lattices are extreme 
narrow-band systems: for instance, in the {L, i? = 5} 
family the conduction-electron band width Wc is in the 
range 0.11 - 0.14 eV; in the {L,R = 7} family Wc is 
even smaller (0.02 — 0.04 eV). The band gaps, unlike band 
widths, decrease with L for given R. In the {L,R ~ 5} 
family, they are in the range between Eg — 0.146 eV 
(for L = 19) and = 0.735 eV (for L = 9); in the 
{L, R ^ 7} family, the corresponding range is between 
Eg ^ 0.118 eV (for L = 20) and Eg = 0.301 eV (for 
L — 12). The band gap sca ling with antidot dimensions 
was studied in Refs. [ij and[lj. 

It is worth noting that as a consequence of the bipartite 
nature of the underlying honeycomb lattice, the resulting 
energy spectrum of our nearest-neighbor tight-binding 
model displays particle-hole symmetryjii This property 
is not perfectly retained in the exact band structureJ^ 

IV. PHONON SPECTRA OF GRAPHENE 
ANTIDOT LATTICES 

The phonon spectrum of graphene was studied exten- 
sively, using either DFT methods or effective models<^"— 
The effective models are either based on approximating 
directly the interatomic force constants, with terms that 



describe coupling of atoms up to some maximum dis- 
tance, or on adopting an analytic expression for the in- 
teraction energy of two or more C atoms. In either case, 
they contain adjustable parameters whose values are de- 
termined by fitting available experimental- and/or DFT 
data. A well-known class of empirical interatomic po- 
tentials for carbon are the Tersoff-Brenner potentials.— 
However, these potentials, with interatomic-interaction 
range that corresponds only to second-nearest-neighbor 
distance, do not provide the accuracy typically needed 
in phonon calculations i^i*^ A generalization of these po- 
tentials, with interactions up to fourth- nearest-neighbor, 
was implemented by Tewary and Yang for calculating 
the phonon spectra and elastic constants of graphene<^ 
Yet, this graphene-specific generalization does not lend 
itself to the use in other carbon-based systems, as the 
generalized potentials become unstable away from the 
equilibrium lattice configuration of graphene. 

In the present work, we determine the phonon spec- 
tra of the {L,R — 5} and {L,R = 7} graphene antidot 
lattices (with the values of L specified above) using two 
independent models that have recently been shown to 
yield very accurate results for graphene itself: the fourth- 
nearest-neighbor force-constant (4NNFC) method^^ in 
the parametrization of Zimmermann et al.^ and the va- 
lence force-field (VFF) model, developed for graphene by 
Perebeinos and Tersoff^^I These two models belong to the 
two aforementioned mutually distinct groups of effective 
lattice-dynamical models. Both models are, in principle, 
applicable to any sp^-hybridized carbon system. 

The 4NNFC model, introduced by Saito et alM, en- 
tails a direct parametrization of the real-space force con- 
stants up to fourth-order neighbors, hence the name. The 
model includes a set of 12 adjustable parameters, whose 
values are determined by fitting the ab-initio phonon dis- 
persion (for the parameter values, see Table I in Ref. ,32). 
These force-constant parameters correspond to the radial 
(bond-stretching), and both in-plane and out-of-plane 
tangential (bond-bending) directions between an atom 
and its n-th nearest neighbors. The underlying formal- 
ism is described in detail in Ref. [H. 

The VFF model of Perebeinos and Tersoff is based on 
an explicit expression for the interaction energy, which 
includes six different contributions, respectively related 
to: bond-stretching, bending, out-of-plane vibrations, 
misalignment of neighboring tt orbitals, bond-order, and 
coupling between bond stretching and bond bending [see 
Eqs. (1) and (2) in Ref. [s^. Each contribution is pa- 
rameterized by the corresponding stiffness constant; their 
values (see Table I in Ref. |33|) are determined by fitting 
experimental and ab-initio data for phonon energies and 
elastic constants. The model makes no reference to any 
underlying crystal structure; the only restriction to its 
use is that the local geometry be consistent with sp^ 
bonding, i.e., that the three adjacent C atoms are not 
too far from being 120° aparti^ Thus the model can be 
directly applied not only to graphene but also to various 
other allotropes of carbon, with graphene antidot lattices 
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being an example of such structures. 

In each particular case, we first find the equilibrium 
lattice configuration by relaxing the atoms until the 
forces on all of them are smaller than 10^^ eV/A. We 
then construct the force-constant tensor 



(4) 



where Mm^(R) is the displacement - in direction (3 {{3 — 
x,y,z) from the equilibrium position - of an atom at 
R+dm, and i?tot({u,„(R)}) the total lattice potential en- 
ergy. The normal- mode frequencies WA(q) and eigenvec- 
tors v^(q) = [v^(q), . . . , v^_^^(q)]^ (with A = 1, . . . , STV^t 
enumerating phonon branches) are obtained by solving 
numerically the secular equation for the dynamical ma- 
trix D(q) EE J2r D(R) e-'i-^. 

The salient feature of the computed phonon spectra is 
that the highest optical-phonon energy at q = is essen- 
tially inherited from graphene itself and thus only weakly 
dependent on L and R. Generally speaking, the com- 
paratively high optical phonon energies in carbon-based 
systems (graphene, nanotubes, intramolecular modes in 
fuUerides) - which typically extend up to about 200 meV 
- are related to the low atomic mass of carbon and to 
the stiffness of the C-C bond. The energies we obtain 
for this highest-lying mode in antidot lattices are around 
195.3 meV in the 4NNFC approach and 197.5 meV in the 
VFF approach. The highest-energy phonon modes have 
in common that they do not involve significant displace- 
ments in the vicinity of the antidot edges. 

Given the number of phonon modes involved, the most 
meaningful way of comparing the results of the two meth- 
ods employed entails calculation of the phonon density- 
of-states 



'^A(q)] 



(5) 



q,A 



This quantity is normalized here such that when inte- 
grated over the whole range of cj it yields the total 
number of the degrees of freedom per unit cell, i.e., 
the number of phonon branches 3iVat. In particular. 
Figs, [^fa) and ^h) illustrate such a comparison for the 
{L = 9, R = 5} and {L = 13, i? = 7} antidot lattices, 
respectively. The plots show the qualitative agreement 
between the two approaches, as well as very similar char- 
acter of the phonon density-of-states in both families of 
antidot lattices studied. The accuracy of both meth- 
ods can potentially be improved as system-specific ex- 
perimental data become available, allowing to determine 
more accurately the values of adjustable parameters. 



V. ELECTRON-PHONON COUPLING IN 
GRAPHENE ANTIDOT LATTICES 

The dominant e-ph coupling mechanism in 7r-electron 
systems is the modulation of hopping integrals due to lat- 
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FIG. 2: The phonon density-of-states for the (a) {L — 9, R = 
5} and (b) {L = IS, R = 7} graphene antidot lattices, ob- 
tained using the 4NNFC (solid line) and VFF (dashed line) 
methods. 



tice distortions,— i.e., Peierls-type coupling.™ The lat- 
ter forms the basis of the Su-Schrieffer-Heeger (SSH) 
modeli^i^Z In sp^-bonded carbon-based systems optical 
phonons modify the length of the in-plane a bond be- 
tween two C atoms, thus altering the overlap between 
the out-of-plane tt orbitals centered on these atoms. As a 
result, the 7r-electron hopping integrals are dynamically 
bondlength-dependent; to a linear approximation, they 
are proportional to the bondlength modulation, with 
a = 5.27 eV/A being the corresponding e-ph coupling 
constanti22i2S The coupling to acoustic phonons, on the 
other hand, is known to be rather weak in graphene 
and carbon nanotubes i^^'^^ In our system, where acous- 
tic phonons have even lower energies than in graphene 
and carbon nanotubes, this coupling is expected to be 
even weaker and can therefore be safely neglected. 

The momentum-space form of the e-ph-coupling 
Hamiltonian reads 



1 



E 



{bl X + &q,A) 



k,q,A,?T,n' 



(6) 

where a„_k (a|j ^J destroys (creates) an electron with 
quasimomentum k in the n-th Bloch band, and 7^„, (k, q) 
stands for the (bare) e-ph interaction vertex function. 
Quite generally, j^^, (k, q) has the meaning of the ma- 
trix element - between the electronic Bloch states (n', k) 
and (n, k -f q) - of the induced self-consistent potential 
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AL^Aq per unit displacement along the phonon normal 
coordinate that corresponds to the mode Aq: 



iL'iK q) = 



1 



i',k|AC/Aq|n,k + q) . (7) 



In systems where AU\q can be determined as a by- 
product of a DFT band-structure calculation (in which 
case it represents a derivative of the self-consistent Kohn- 
Sham potential), like for graphenC) ^^'^^ the last expres- 
sion incorporates contributions of the phonon mode Aq to 
all possible e-ph coupling mechanisms. Hence there is no 
need to assume a particular microscopic e-ph interaction 
form. Here, however, where such an approach is unfea- 
sible, evaluation of 7,^„, (k, q) is based on the adopted 
real-space form of the e-ph coupling written in the tight- 
binding electron basis. In the case at hand 



7l,(k,q) = l/^„,(k,q) + VFl,(k,q) 



(8) 



where V^^, (k, q) and W^„, (k, q) are the respective con- 
tributions of the neighbors within a single unit cell and 
those from adjacent unit cells. As shown explicitly in 
Appendix 



^e„,(k,q) ^ [vi+,(q) - vi(q)] 

X + , (9) 

where 6 = 5/||5||, while m + S denotes neighbors d,,,, + S 
of site d„i, and 



W^l'(k,q) = 



= S- [e'^-^i, (q)- v,^,(q)] 

r _ — ?(k+q) ■a//^n.k+q \ * /^n' .k i ^-ik-a //^n.k+q \ * y^n .kl 

X [e K^mi ) ^m' 'T-e [ly^ ) J . 

(10) 

The prime on the last sum indicates that the sum is re- 
stricted to the neighbors d„i + S oi site that satisfy 
the condition d,,,, + ^ — a + for some mi = mi{d), 
with a = ±ai, ±a2, ±(ai — a2) (see Fig. [IJ. 

For the conduction band (n — > c) electrons, in the ab- 
sence of interband (n ^ n') scattering, the e-ph Hamil- 
tonian of Eq. ([5]) reduces to the effective form 

^ip =^T. 7c'c(k,q)aU+qae,k(6^_q,, + Sq,A). (11) 

k,q,A 



Evaluation of the e-ph vertex functions 7^(.(k, q) based 
on Eqs. ([9]) and ([10]) requires the full knowledge of the 
phonon polarization vectors v^^ (q) and Bloch wave func- 
tions, with the latter ones entering through the coeffi- 
cients C2='='^. 



VI. PHONON-INDUCED RENORMALIZATION 

A. Consequences of the momentum dependence of 
e-ph coupHng 

At variance with most other familiar types of e- 
ph interaction, which are either completely momen- 
tum independent (Holstein-type coupling22i) or only 
weakly phonon-momentum dependent (Frohlich-type 
coupling^), Peierls-type coupling depends strongly on 
both the electron and phonon momenta. This translates 
into a strongly-retarded nature of the interaction, which 
can have a variety of physical consequences. ^-'^ Due to 
particular geometry of the system, the momentum depen- 
dence of e-ph interaction in the case at hand is clearly 
more complicated than that of the standard SSH cou- 
pling, which is usually studied on monoatomic lattices 
(one-dimensional, or a two-dimensional square one). 

The quantity that characterizes phonon-induced renor- 
malization for conduction-band electrons is the quasi- 
particle spectral weight due to e-ph coupling, Zc(k) = 



K^ckl-iAck)!^ , where iV'ck) 



|0)e is the bare 



conduction-electron Bloch state (i.e., a common eigen- 
state of He and the total electron momentum operator 
= J2n',k' k'a]^, k'^n'k') ^ud l^ick) that of the coupled 
e-ph system (i.e., an eigenstate of the full Hamiltonian of 
the coupled e-ph system and the total crystal momentum 
operator K = Y.n',k' ^' ai',k'«"'.k' + EA,q' q' ^l,q'^A,q')- 
In a recent work)2S we evaluated Zc{0) for the {L, R — 
5} and {L,R = 7} families of lattices, finding a rather 
strong phonon-induced renormalization Zc{0) — 0.20 — 
0.27 [i.e., Z^^{0) = 3.7 — 5], compared to graphene^^ 
where Z = 0.93 at Dirac points, or largeri^ These re- 
sults, with an excellent agreement between the 4NNFC 
and VFF approaches, seem to suggest that charge car- 
riers acquire polaronic character. In what follows, we 
analyze in detail the most important phonon modes and 
momentum dependence of their corresponding e-ph ver- 
tex functions. 

The largest contribution to the spectral weight (75 — 80 
percent in total) comes from the high-energy modes, 
more precisely two groups of modes centered around 173 
meV and 194 meV (cf. Fig. [5]). These modes have mixed 
character, being neither purely transverse nor longitudi- 
nal. We also note that within the phonon energy intervals 
that contribute to the spectral weight, the contribution 
is dominated by a few modes only. For example, in the 
case of {L = 13,R = 5} [{L = 15,i? = 7}] lattice only 
9 (24) modes are sufficient to provide already 50 percent 
of the spectral weight. 

The highest-energy phonon branches yield the largest 
individual contribution to Z^^{Q) in all the lattices stud- 
ied; these contributions are in the range 0.57 — 0.65 in the 
{L, i? = 5} lattices and 0.27-0.34 in the {L, R=7} lat- 
tices. Examples of the atomic displacement patterns of 
such modes for {L = 13, i? = 5} and {L = 15, i? = 7} an- 
tidot lattices are given in Figs.|3]and|31 from which we can 
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FIG. 3: (Color online) The atomic displacement pattern 
corresponding to the zone-center phonon mode at energy 
Ex ~ 197.5 meV of the {L = 13, J? = 5} graphene antidot 
lattice (top) and its zoomed version (bottom). The displaced 
(undisplaced) atomic positions are designated by the squares 
(dark circles). 



FIG. 4: (Color online) The atomic displacement pattern 
corresponding to the zone-center phonon mode at energy 
Ex f» 197.5 meV of the {L = 15, R — 7} graphene antidot 
lattice (top) and its zoomed version (bottom). The displaced 
(undisplaced) atomic positions are designated by the squares 
(dark circles). 



infer that these modes do not involve significant atomic 
displacements in the vicinity of the antidot edges, but 
only away from them. This is a consequence of the fact 
that the system is "stiffer" to the displacements of atoms 
away from edges (which have three neighbors) than those 
along the edges (only two neighbors). Similar argument 
explains the dominance of tangential displacements (over 
radial ones) in the highest-energy phonon modes of ful- 
lerides."*^ 

The last conclusion about the importance of the 
highest-energy optical phonon branch is similar to that 
in graphene: based on a first-principles study of the e-ph 
interaction in graphene. Park et al.^ concluded that the 
results could be reproduced by assuming a single Einstein 
mode at around 200 meV. 

On the other hand, the optical phonons at energies 
below 30 meV are far less important than their high- 
energy counterparts, with their cumulative contribution 
to the overall spectral weight being at most 20 percent 
in the antidot lattices studied. This can already be un- 
derstood based on the shape of the phonon density-of- 
states (Fig. [2]) , from which it is obvious that in the low- 
energy region the phonon density-of-states is compara- 
tively small. 

In order to get a qualitative insight into the obtained 
results, it is instructive to analyze in detail the underly- 
ing e-ph coupling mechanism. The momentum depen- 
dence of the bare vertex functions for an electron at 
the conduction-band bottom (k = 0) and the dominant 



phonon modes in typical lattices is depicted in Figs.[51Ja)- 
(d). What can be inferred from the plots is that the 
e-ph coupling has a very rich momentum dependence, 
characterized by a high degree of anisotropy within the 
irreducible wedge of the BZ. We find several characteris- 
tic momentum-space profiles for the coupling functions, 
some of which correspond to e-ph coupling peaked at (or 
around) q = [Figs. [SIa),(c) and (d)], and the remain- 
ing ones to the cases where the latter have maxima at 
nonzero phonon momenta [Fig. [5l^b)]. 

In all the cases considered, the coupling to the highest- 
energy phonon branch is strongest at zero phonon mo- 
mentum and in its immediate vicinity. The last ob- 
servation leads us to conclude that the small-phonon- 
momentum (q « 0) scattering involving this mode is 
to a large extent responsible for the sizeable phonon- 
induced renormalization that we have obtained. As can 
be inferred by analyzing the expression for Z~^{0) in 
the lowest-order Rayleigh- Schrodinger perturbation the- 
ory [Eq. (8) in Ref. [23, the q w terms yield the largest 
contribution to the relevant integral since for small q val- 
ues the denominator of the integrand is very small. The 
last assertion about the relevance of q « scattering is 
quite general and holds even when the e-ph interaction is 
completely momentum independent (local in real space), 
as is the case for Holstein-type interaction. However, our 
system constitutes an extreme realization of such a phys- 
ical scenario as the function |7cc(k = 0,q)|^, correspond- 
ing to the highest-lying phonon, attains its largest values 



precisely at q = and thus leads to a large mass enhance- 
ment. Needless to say, distinct examples of momentum- 
dependent e-ph interactions do exist where the physical 
situation is rather different; for the SSH coupling^iiil 
7ssH(k, q) cx sin(k • a) — sin [(k + q) • a] one has 7ssH(k = 
0, q) oc |q| at small q. Thus the SSH coupling for k = 
vanishes as q — 0. Another familiar example is furnished 
by the coupling of Zhang- Rice singlets to the "breathing" 
(oxygen bond-stretching) phonon modes in cuprates^ 

where 7(k, q) = 7(q) oc sin^(ga;/2) -|- svo?{qy/2) and is 
therefore also suppressed at small phonon momenta. 




FIG. 6: Pictorial representation of the Dyson equation (top). 
lUustration of the SCBA self-energy, given by a sum over all 
noncrossing diagrams (bottom). 



B. Nonzero- momentum spectral properties due to 
e-ph coupling 

In the following, we study nonzero-momentum spectral 
properties due to e-ph coupling. Because of a consider- 
able computational burden involved, we perform calcu- 
lations on selected antidot lattices - {L = 9,i? = 5} and 
{L = 13, i? = 5} from the {L,R = 5} family, as well as 
{L = 15, i? = 7} from the {L, R = 7} family - within the 
VFF approach. 

We determine the conduction-band electron self-energy 
using the self-consistent Born approximation (SCBA) 4^ 
In the SCBA the electron self-energy is obtained by sum- 
ming over all the non-crossing diagrams, as illustrated in 
Fig.[6l Although this approximation neglects e-ph vertex 
corrections in the self-energy calculation, it yields good 
results for not-too-strong e-ph coupling4i The SCBA 
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FIG. 5: (Color online) The q-dependence of the moduli 
|74(k = 0,q)| of the bare e-ph vertex functions, for selected 
high-energy phonon branches A, in the hexagonal Brillouin 
zones of the {L = 13,R = 5} [(a)-(b)] and {L = 15,7? = 7} 
[(c)-(d)] graphene antidot lattices. The energies of the cor- 
responding optical phonons at q = 0, as obtained using 
the VFF approach, are Ex « 197.5 meV [(a),(c),(d)] and 
Ex « 197.4 meV [(b)]. 



self-energy is given by 



Sc(k,-)=2^E/^"l7cc(k,q: 



xDxiq,n)Gc(k + q,LU-n), (12) 



where 



Gc(k,L.) 



1 



(13) 



uj — Ec (k) — Sc (k, w) + iS 

D^{q,n) ^ ^— r^- ^^ ^ (14) 

it — ujx + iO+ il + LUX — iO 

are the interacting electron and the free phonon propaga- 
tors, respectively. The frequency integration in Eq. (jl2p 
can be carried out explicitly, leading to 



^c(k,^)-^E l7cc(k,q)|^ 



Lj - LJx- £c(k -I- q) - Sc(k -I- q, w - ux) + iS 



(15) 



The last equation is solved self-consistently for elec- 
tronic (conduction-band) self-energy I]c(k, oj) using the 
method of iteration, with the self-energy from Rayleigh- 
Schrodinger perturbation theory^S taken as an initial 
guess in the iterative process at every k and oj. The 
resulting SCBA self-energies for the {L = 13, R = 5} 
and {L — 15, R — 7} antidot lattices are depicted in 
Fig. Elja) and (b), respectively. In this self-energy calcu- 
lation, a broadening parameter 6 was used, with a value 
of 5 meV that corresponds to the resolution of state-of- 
the-art photoemission experiments.— The SCBA elec- 
tron Green's function is then evaluated based on Eq. , 
with a typical example shown in Fig. [51 The obtained 
results are to be compared with future ARPES measure- 
ments of the single-particle spectral function Ac(k, a;) — 
— ImG'c(k, ll')/tt. 

In the case of very weak e-ph interaction the peaks in 
the spectral function appear at energies that are shifted 
from the bottom of the dressed-electron band by a char- 
acteristic phonon energy. In the case at hand, the po- 
sitions of peaks result from a complex interplay of mul- 
tiple phonon modes involved and momentum-dependent 
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e-ph coupling. The fact that the second peak is shifted 
from the first one by an energy larger than the maximal 
phonon energy in the system is an additional signature 
of strong e-ph interaction and polaronic behavior. 

The renormalized conduction-band electron dispersion 
i?c(k) is obtained by solving self-consistcntly the equa- 
tioni^ 



£;c(k) = £c(k) -f ReEc[k,i;c(k)] , 



(16) 



where ReEc(k, w) is the real part of the SCBA self- 
energy. The quasiparticle weight .^c(k) is then computed 
for different nonzero quasimomenta k using the standard 



expressioi 



— ReE,(k,w)|^^^^(^^ 



(17) 



The obtained values for the {L — 9, R — 5} antidot lat- 
tice, for example, are in the range between 0.33 and 0.37, 
all of them being larger than Zc(k = 0) w 0.208 ob- 
tained using Rayleigh-Schrodinger perturbation theory. 
The increase of the spectral weight away from k = 
is intimately related to the fact that at larger total po- 
laron momentum the momenta of phonons relevant in 
the scattering process are also larger. Accordingly, the 
e-ph coupling that predominantly decreases with phonon 
momentum (recall Sec. IVI Al) leads to a smaller renor- 
malization. 
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FIG. 8: (Color online) The SCBA electron propagators in 
the (a) {L = 13, R = 5} and (b) {L = 15, i? = 7} graphene 
antidot lattices for k — 0.2 b, where b = (0, 1) * 47r/(3I/a). 
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FIG. 7: (Color online) The real and imaginary parts of the 
SCBA conduction-band electron self-energy in the (a) {L — 
13, -R = 5} and (b) {L = 15, i? = 7} graphene antidot lattices 
for k = 0.2 b, where b = (0, 1) * 4n/{3La). 



In order to characterize the phonon-induced renor- 
malization, it is of interest to determine the direction- 
dependent e-ph mass enhancement factor 



1 {a = x,y), 



(18) 



where rrios is the effective (in the presence of phonons) 
and m* the bare-band mass. Starting from the defin- 
ing self-consistent relation for the renormalized disper- 
sion [Eq. ([T6|) ]. it is straightforward to show that (see 
Appendix [B] for a derivation consistent with our present 
notation; alternatively, see Ref. |47| ) 



d 



1-— ReEe(k,..)|,^„__^^ 



(0) 



decCka, 



■ ReEc(kQ, w)| 



(19) 

where k^ = (k • ea)ea = kaSa- It then directly follows 
that 



^c-Ho) 



d 



dedka) 



Rei;c(ka,a;) 



-1. (20) 



The last quantity can be obtained by a numerical dif- 
ferentiation of ReEc, using previously obtained values 
for Z~^{0). For the {L = 9, i? = 5} lattice, our calcu- 



lations yield the values A 



2.411 and A 



(y) 



2.448 
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[the calculated values of the respective denominators of 
the fraction in Eq. ^ are 1.4106 and 1.3954]. For the 

{L ~ I'S, R ~ 5} lattice we obtain Amc = 2.389 and 
X^^l = 2.081, respectively. In the {L = 15, i? = 7} lattice 
case, we get Ami = 1.365 and Ame = 1.255. 

The obtained numbers suggest that the anisotropy in 
the e-ph mass enhancement is relatively weakly pro- 
nounced, despite the existence of anisotropy in the 
momentum-dependent e-ph coupling [recall Fig. [5] . This 
can be qualitatively explained by observing that the lat- 
ter anisotropy of the e-ph coupling is smallest for the 
highest-energy phonon modes which yield the strongest 
coupling [for an illustration, note the variation of the val- 
ues of |7cc(k = 0,q)| in Fig. [Sfa) and (c), as opposed 
to those in (b) and (d)]. We can therefore draw the 
conclusion that the overall effective-mass anisotropy in 
graphene antidot lattices is predominantly determined 
by the anisotropy of the bare band mass itself, rather 
than that stemming from phonon-related effects. 



VII. DISCUSSION 

The strong renormalization obtained suggests that the 
charge carriers in the system acquire polaronic character. 
Indeed, it is not surprising to have polaronic charge car- 
riers in a narrow-band system. This is a commonplace, 
for example, in organic semiconductors, such as organic 
molecular crystals, where narrow electronic bands result 
from weak van der Waals intermolecular bonding and the 
e-ph interaction has similar (slightly nonlocal, but still 
short-ranged) characterJ^ Compared to these systems, 
graphene antidot lattices have yet narrower conduction 
bands (e.g., bare conduction-electron band widths in the 
polyacene organic crystals can be as large as 700 meV) 
and lower dimensionality, both of these factors being con- 
ducive to the larger mass enhancement found here. 

The strong-coupling regimes in interacting e-ph sys- 
tems are usually explored using exact diagonalization 
or variational methods. The criteria for the onset of 
a (smooth) polaron crossover with the increasing e-ph 
coupling strength are typically based either on the be- 
havior of the average number of phonons in the polaron 
ground state - which in the crossover region grows precip- 
itously - or on the behavior of the quasiparticle spectral 
weight, which in the same regime undergoes a sharp de- 
crease. For a system with short-range, non-polar e-ph 
interaction (regardless of its concrete form), as is gener- 
ally the case in covalently-bonded materials such as the 
one considered here, the polaron crossover has the nature 
of a change from a quasi-free electron to a small polaron 
(rather than being a large-to-small polaron crossover, for 
which a long-range character of the e-ph interaction is a 
crucial ingredient). However, for a system of such com- 
plexity as studied here the use of the aforementioned 
methods would clearly be inconceivable even with a sin- 
gle "effective" phonon mode, thus the criteria for pola- 



ronic behavior cannot be utilized in a strict sense. The 
perturbative approach is the most powerful technique at 
our disposal. Obviously, the obtained results are not ex- 
pected to be valid quantitatively. However, they should 
still be qualitatively correct, providing an unambiguous 
evidence for the relevance of phonons in graphene antidot 
lattices. In addition, it should be stressed that the extent 
to which the obtained results can be improved is limited: 
while the accuracy in obtaining the phonon spectra can 
certainly be improved, due to the prohibitive system size 
the band structure remains tractable only at the level 
of the tight-binding method. Finally, even if DFT cal- 
culations on such a complex system were feasible, the 
validity of DFT when calculating e-ph coupling in two- 
dimensional systems has been called into question4^ 

An experimental verification of the polaronic nature of 
carriers in our system should be possible in future opti- 
cal absorption measurements. Some qualitative predic- 
tions can already be made based on general features of 
(intraband) optical absorption in systems with Peierls- 
type (SSH) e-ph coupling. In the polaronic regime, the 
ground state of such systems is characterized by a short- 
ened bond - being a result of two adjacent sites shifting 
towards one another - on which the electron is local- 
ized. The ground state is then the even combination 
of the two local electron states. Viewed through sim- 
ple Franck-Condon-type picture, the optical absorption 
in such systems proceeds through two different channels: 
one in which the electron is excited onto another bond 
(which is not shortened), and the other one in which 
the electron is excited from the (even-symmetry) ground 
state to the local odd-symmetry state on the shortened 
bond itself. The excitation energy in the first process 
corresponds to the energy gained by the lattice distortion 
(compared to the undistorted one) and is given by 2i?p, 
where Ep is the polaron binding energy. In the second 
absorption channel, which corresponds to a transition be- 
tween a state where the electron energy is lower by 2Ep 
with respect to the undistorted lattice to a state where 
it is higher by the same amount, the excitation energy is 
twice as large (4i?p). When the finiteness of the relevant 
phonon frequency is taken into account, the absorption 
peaks at energies 2Ep and 4Ep from the Franck-Condon 
picture broaden into bands characterized by phonon fea- 
tures separated by the characteristic phonon energy. 

An immediate implication of our results is that the 
charge transport in graphene antidot lattices - at vari- 
ance with that of their "parent material" - ought to be 
modelled by taking into account the inelastic degrees of 
freedom. This was done previously with carbon nan- 
otubes,^^ where the inelastic effects in transport are rele- 
vant only in the high-bias regime as a consequence of the 
carrier energy exceeding that of an optical phonon (with 
energy of « 200meV)fi^ However, in our system there 
are also optical phonons at far lower energies. Thus the 
inelastic effects should be even more pronounced than in 
carbon nanotubes and not necessarily restricted to the 
high-bias regime. 
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It is worthwhile noting that the recently fabricated 
graphene antidot lattices*' comprise even larger num- 
ber of atoms per unit cell (iVat ^ 5,000) than the ones 
discussed here. For comparison, their periods are in the 
range 27—39 nm, while the largest period of the superlat- 
tices studied here is around 8.2 nm. Thus it is clearly of 
interest to extend the present study to such cases. How- 
ever, a detailed treatment of the e-ph coupling effects, 
whose complexity scales roughly with N^^. , for systems of 
such size would be extremely computationally demand- 
ing even with the methodology employed in the present 
work. Such cases can, however, be addressed within the 
effective-mass approximation for electronic band struc- 
ture (wherein the wave functions are expressed through 
slowly-varying envelope functions) and with a contin- 
uum treatment of phonons (deformation potential), an 
approach well established in other systems.— Knowing 
that these large-period lattices have larger conduction- 
band widths (in addition to having smaller band gaps) it 
is likely that the e-ph coupling effects will be quantita- 
tively less important than in the lattices studied here. 

Another important issue for future investigation is 
the interplay of (short-range) Peierls-type coupling in 
graphene antidot lattices with the long-range (Frohlich- 
type) e-ph coupling at their interface with polarizable 
substrates such as SiC or Si02.^^ This interplay between 
the two e-ph coupling mechanisms is likely to be relevant 
in a field-effect-transistor geometry, where a graphene 
antidot lattice can be used as the semiconducting chan- 
nel and Si02 as a high-K gate dielectric* Quite recently, 
motivated by its relevance in organic field-effect transis- 
tors, such interplay has been investigated within a sim- 
plified one-dimensional model. The main conclusion of 
this study was that even a rather weak SSH coupling is 
sufficient to obtain a polaronic behavior in the simultane- 
ous presence of a moderate polar coupling. In graphene 
antidot lattices, as a consequence of their different ge- 
ometry and the ensuing momentum dependence of the 
e-ph coupling (compared to the genuine SSH coupling; 
recall discussion in Sec. IVI A|) such issues would need to 
be reconsidered. 



VIII. SUMMARY AND CONCLUSIONS 

To summarize, we have studied the effects of electron- 
phonon interaction in graphene antidot lattices. Using 
realistic band-structure and phonon-spectra calculations 
as an input, we have described the system based on a 
model that accounts for the phonon-modulation of elec- 
tronic hopping integrals - Peierls' type electron-phonon 
coupling. We have demonstrated a significant electron- 
phonon mass enhancement in this system, which is likely 
to be verified in future ARPES measurements. 

The main message of the present work is that in 
graphene antidot lattices - in contrast to graphene, 
their parent material - optical phonons, particularly the 
highest-energy ones, play a very important role. This can 



be attributed to the narrow-band character of these sys- 
tems, and, especially, to the peculiar momentum depen- 
dence of the relevant electron-phonon interaction which 
is strongest for small phonon momenta. 

Apart from being related to a particular family of 
graphene-based nanostructured systems of keen current 
interest our study also bears some fundamental im- 
portance. Namely, it provides an insight into the rele- 
vance and physical consequences of strongly momentum- 
dependent electron-phonon interactions. The body of 
work concerned with such interactions is fairly limited, 
even at the level of model investigations, in spite of 
their importance in realistic electronic materials. 

The present work bodes well for further studies of 
graphene antidot lattices, in the finite-temperature case 
and in the presence of a substrate and/or metallic con- 
tacts,— all of these aspects being of particular relevance 
for future room-temperature, carbon-based electronic de- 
vices. By providing evidence for the relevance of the 
inelastic degrees of freedom, our findings are expected 
to also facilitate future studies of electronic transport 
in these systems. Detailed investigations of electron- 
phonon-interaction effects in other emerging classes of 
graphene-based superlattices^ are clearly called for. 
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A. DERIVATION OF THE 
ELECTRON-PHONON-COUPLING VERTEX 
FUNCTIONS 

By introducing the matrix Ck (for each k in the BZ) 
such that C,"''' = (Ck)njn, Eq. ^ can be rewritten as 

tAk(r) = Ck0k(r) , (Al) 

where i/'k(i") = [V'nk(r)| n = 1, . . . , N^.t]'^ and (f>^{r) = 
[0mk(r)| m — 1, . . . , Ni.t]'^ are A^at-component vectors. 
The inverse of the last equation 

0k(r) = C^i^k(r) , (A2) 
when recast componentwise reads 

0mk(r) = ^(Cki)„„i^„k(r) . (A3) 

n 

Given that C,^'*' are the coefficients of a linear transfor- 
mation [recall Eq. ^] between two complete orthonor- 
mal sets of functions, the matrix Ck is unitary, so that 
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k }ran 



is equivalent to 



{Cl^^^y. Consequently, Eq. (|M 



,k(r) = 5](C;;-'')*V^«k(r) 



(A4) 



and, combined with the inverse Fourier transformation 
^(r-R-d„O = ^V0™k(r)e-*-^, (A5) 

readily implies that 

4+d„ = ;4 E(^™')*^^^' ''4k ■ (A6) 

n,k 

If, in addition, we assume that for a given vector d,„, 
the vector d,„ + S denotes a site within the same unit 
cell (i.e., hopping takes place between two lattice sites 
within the same unit cell), it also holds that 



n,k 



(A7) 



From the last two equations, combined with the expres- 
sion for phonon displacements 



UA,R+d„ 



<{<l), (A8) 



we readily obtain 

y^(°R.+d„,+a"R+d„i +H.C.) [uA.R+d„+g -UA,R+d„J 



R 



A^„,^k,q V2M^A(q) 



^•[vi+,(q)-v^(q)] 



r/,^n,k+q\*^n',k , f /^n,k+ci\* ^n' M 1 

q.A + ^q,A) • (A9) 

The last expression straightforwardly leads to the part 
of the Hamiltonian in Eq. ^ that corresponds to the 
hopping within the same unit cell, from which we can 
read out the expression for V^^, (k, q) in Eq. (0) . 

If, on the other hand, for a vector in the given unit 
cell, the nearest-neighbor vector dm -I- S belongs to one 
of the adjacent unit cells - displaced from the given one 
by the vector a [one of the vectors ±ai, ±a2, ±(ai —32)] 
- Eq. (|A7|) ought to be replaced by 



aR+d„+5 = 2^ (-'mi e 

^ n.k 



k-(R+a) 



'^nk 



(AlO) 



where in the last expression mi = mi{S) (mi e 
{1, . . . ,iVat}) is chosen in such a way as to satisfy the 



condition d™ -|- 5 = a -I- dm^ . In an analogous manner, 
the expression for UA,R+d,„+<5 ~ UA,R+a+d„, should then 
read [cf. Eq. (jA8|) ] 



jq-(R+a)/Lt 



UA,R+d,„+5 



N ^ 
q 



- (q) 



v/2AfcjA(q) 

(All) 

Now, by repeating the steps in the above derivation of the 
expression for (k, q), with regard for Eqs. (|A10[) and 
()All[) . we easily recover expression (fTOl) for W^„/(k, q). 



B. DERIVATION OF EQ. ([T9|l 

We start from the defining self-consistent relation^^ for 
the renormalized conduction-band electron dispersion: 



£;e(k) =£c(k)-hReI]e[k,i;,(k)] 



(Bl) 



By expanding the second term on the right-hand side of 
the last equation around k = 0, we readily obtain 



d 



(k) = (k) + k . — Re Ee (k, .;) | ^^^^^^^^ 



(0) 



d 



i?e(k)^ReI]e(k,a.)|,^„_^^(„) +0(k2) , (B2) 



which combined with the identity 

Se(0) -ec(0) + ReEc[0,£;c(0)] 
leads to the relation 



(B3) 



[£;c(k) -£;c(o)] 



1- a^i^eEc(k,^)|^^P^^^^^(g^ 



= ec(k) - ec(0) + k • ^ ReEc(k, w)|i^^o_^^^^(o) 

-h Oi^) . (B4) 



The mass renormalization in direction a is given by 

£c(ka)-ec(0) 



TOcff 



— lim 



(B5) 



where k^ = (k • ea)ea = kaGa. Using the special case 
k = ka of Eq. (jB4[) . after some straightforward manipu- 
lations, we finally obtain 



d 



m* 



l-ajR°Sc(k,a;)|^^„^ 



i5c(0) 



d 

^^-^ReEc(k„,a;)|^^^Q_^^^^(o) 

(B6) 
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